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TECHNICAL NOTE 2165 

METHOD OF ANALYSIS FOR COMPRESSIBLE FLOW THROUGH MIXED -FLOW 

CENTRIFUGAL IMPELLERS OF ARBITRARY DESIGN 

By. Joseph T. Hamrick, Ambrose Ginsburg 
and Walter M* Osborn 


SUMMARY 

A method is presented for analysis of the compressible flow 
between the hub and the shroud of mixed -flow impellers of arbitrary 
design* Axial symmetry was assumed, but the forces in the meridional 
(hub to shroud) plane, which are derived from tangential pressure 
gradients, were taken into account. 

. The method was applied to an experimental mixed -flow impeller. 
The analysis of the flow in the meridional plane of the impeller 
showed that the rotational forces, the blade curvature, and the 
hub -shroud profile can Introduce severe velocity gradients along the 
hub and shroud surfaces. Choked flow at the impeller inlet as deter- 
mined by the analysis was verified by experimental results. 


INTRODUCTION 

Because of the complexity of the numerical solution for a 
three -dimensional analysis of flow through a centrifugal compressor, 
most treatments of the problem have been concerned with a two- 
dimensional analysis of the flow between the pressure side of one 
impeller blade and the suction side of the preceding blade. Relax- 
ation methods, in which it was assumed that no velocity variation 
existed between the hub and the shroud, have been used (reference 1) 
to give an analysis of the flow between' blades. La the conventional 
mixed-flow centrifugal impeller, however, the nonsymmetrical hub- 
shroud profile, the nonuniform energy addition across the passage, 
and the blade -angle variation due to close adherence to radial blade 
elements may introduce velocity gradients between the hub and the 
shroud comparable to those between the blades. For mixed -flow 
compressors, therefore, an analysis of the flow in the meridional 
plane is required. (The meridional plane is one that passes through 
and is parallel to the axis of rotation.) 
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Tn order to Investigate the isentropic compressible flow in the 
meridional plane , a flow equation based on equilibrium, cont inui ty, 
anH general energy considerations was derived at the HACA. Lewis 
laboratory and is presented herein* A-rial. symmetry was assumed; but 
the forces that result from tangential pressure gradients were taken 
into account* The numerical solution involves the adaptation of 
stream-filament theory and numerical integration of all equations. 

t Pha applicability of the method for predicting the flow wi thin 
the impeller was investigated by analyzing the flow through a mixed- 
flow centrifugal impeller for which experimental data are available 
(reference 2). Flow streamlines; pressure distributions; and 
velocity distributions relative to the impeller far a given set of 
operating conditions are presented. Theoretical shroud pressures 
and theoretical ma-rimm weight flows over a range of speeds are com- 
pared with those experimentally obtained. 


ANALYSIS 

An analysis is made in the meridional plane of the mixed-flow 
impeller. Axial symmetry is assumed; but forces in the meridional 
pinna that are derived from tangential pressure gradients are taken 
into account. By the use of the methods presented in this a n alysis; 
the isentropic compressible flow in the meridi ona l pla ne may be 
obtained nnd shown in terms of the flow streamlines and the velocity 
and pressure distributions. 

Geometrical quantities . - A complete geometrical description of 
the blade is necessary for a numerical solution. Geometrical quantities 
are defined in appendix A or shown in figures 1 and 2. The blade 
angle 0 is oriented in the following manner: In figure 1(a), a 
line OB, which intersects the axis of rotation at 0, is drawn 
tangent to the point P on a streamline in the maridonal plane. If 
the tangent line OB is rotated about the center line of impeller 
rotation 0-0’ , a conic surface of revolution that intersects the 
impeller blades is formed. This intersection is shown on a developed 
surface of the cone .in figure 1(b). At the point of tangency to the 
streamline (point P), the blade mean line forms the angle -0 with 
the tangent line. . 

Limitations assumptions* - The fluid is assumed to be non- 
viscous but compressible. All compression is by the adiabatic, 
isentropic, steady-flow process. The velocity is assumed uniform 
from blade to bla de in the direction of rotation (axial symmetry) and 
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the average angle of flow relative to the impeller is assumed the 
same as the "blade angle p. This assumption Is considered valid 
where channel flow prevails, hut is questionable at the inlet and 
the outlet* 

Method of analysis . - The method is based an one proposed by 
Flugel for analysis of flow in the blade -to -blade direction and dis- 
cussed by Stodola in reference 3. The forces acting on a particle 
along the streamline path are considered. Hie component of the blade 
surface force normal to the meridional streamline, the centrifugal 
forces resulting from the curvilinear path of the particle, and the 
force due to rotation about the axis of the impeller are equated to 
the change in pressure across the particle in the hub -to -shroud 
direction. The resulting force equation is combined with one based 
an the general energy equation to give an equation involving the 
relative velocity q and a number of variables dependent upon the 
geometry of the impeller. The variables that lend themselves to 
graphical treatment are lumped into two parameters giving a first- 
order linear differential equation linear in q. This equation is 
used in conjunction with the continuity equation to obtain the 
velocity distribution. Streamlines between the hub and the shroud 
similar to those shown in figure 1(a) are drawn in by estimate. By 
solving the differential equation for q and assuming a velocity 
at the hub, the velocity distribution from hub to shroud may be 
obtained by numerical integration. By means of the continuity equation, 
the weight flow corresponding to the assumed hub velocity is obtained. 
Several solutions nay be necessary for obtaining the desired weight 
flow. The streamlines may be rearranged to match the new values of 
velocity and density and the complete analysis repeated. 

Force equation . - For the farces acting on the fluid particle in 
the n direction on the meridional plane (fig. 1(a)), the following 
equation is obtained (appendix B): 

| £ frr + . u ) 2 ar . (1) 

pdnr c rdnn v/ 

From figure 1 it can be seen that 

v = q cos 6 


and 


u = q sin 0 
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In. terms of q., 'equation (1) becomes 


f % = * 2 (- * =£* g)-' i ♦ ** 8ta S I - 

(2) 

TPnaygy addition * - From general energy considerations , the 
isentropic adiabatic energy addition is 

| j^cor! ( cor! + u x ) - toX^J » h x + ^ jour^ + u^ 2 + y x 2 J - % 

where the prerotation term X = r i ( oo r^ + u^ ) 


«2r 2 


1 . V 


h l “ 2g * 2g • 2g 


H i " 


wX-j 

T 


(3) 


From, figure 1, 


g.2 s x? + v 2 


with the assumption that is constant for all streamlines, 

\ - 4 E“ r i )2 • *j*] + ■ “o - ^ W ' -] + 

\ " >0 ■ k f 2 ("l 2 - 'o') ’ («1* - «0*) ' 2 “ ( S. - X o] 


For an isentropic process 


/* 


h - h = 
1 0 


dp 

PS 


Therefore 
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k dp " k ( r i 2 " r o 2 ) - (*i* ' *o 2 ) - 2w 


( x i ■ g] 

(4) 


When differentiated with respect to n, equation (4) hecomes 


ld£ 
p dn 


dr dq d\ 

“ 2r e - « ^ - “ e 


(5) 


General velocity equation . - By combi n i n g equations (2) and (5), 
the following equation is' obtained : 


dn 


t 


cos 2 P 


sin 


_P dr 
dn 


)• 




( 6 ) 


For straight blades lying on the meridional plane with no pre- 


rotation of the fluid at the inlet, the terms 


“n 


and 


0) d\ , 


out; a first-order differential equation linear in q thus remains* 
In order to arrive at a solution for curved -blade Impellers and far 
prerotation at the inlet, it may be desirable to solve far q and 
neglect these terms in the first solution. The values of q thus 
found may be used in the force and prerotation terms far the second 
solution and may be treated in the same manner as those terms 
dependent upon the geometry of the impeller. 

Telocity equation far impeller of arbitrary blade shape having 
radial blade elements with no prerotation . - The blade -farce com- 
ponent F Q normal to the meridional streamlines is evaluated in 

appendix C: 

~ - sin P tan a sin a + 1 sin 2 p sin a tan a 


The velocity equation for a radial -bladed impeller without prerotation 
becomes 


& 

dn 


COS 2 P 


_P dr 
dn 


sin 2 p 


sin a tan a 


q ( : 

2 a) sin P ~ - sin P tan a ^2u> sin a + ) 


)- 


(7) 
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Variables dependent upon the impeller geometry or upon values obtained 
from a previous solution are represented by pare, meters a and b 
vhere 


a 


cos 2 p _ sin 2 P dr _ sin 2 P 

r r dn ” r 

c 


sin a tan a 


and 


b = 2oj sin p ~ + sin P tan a ^2uo sin a + ^ j 
Equation (7) becomes 

When equation (8) is solved for q. 


( 8 ) 



In solving for q, the term du/dm in parameter b is neglected for 
the initial solution in order to obtain estimated values of du/dm for 
the succeeding solution. 


Continuity equation . - The equation for continuity of flow through 
the hub -to -shroud passage may be -written 


w = 2n 


F 

rfpgq cos P dn 

Jo 


( 10 ) 


■where f is a factor that takes into consideration the area taken up 
by the blade cross section. The value of the density P is deter- 
mined from the energy equation (equation (3)) and the 

1 

relation ■ p 

p t,i 
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NUMERICAL EXAMPLE 


Impeller Example 

A numerical example has "been computed for an 18-blade mixed-flow 
centrifugal impeller, the performance of which is described in refer- 
ence 2. The Impeller was designed for a constant cross-sectional flow 
area from inlet to outlet. The blade curvature is such that the 
intersection of the blade with any circular cylinder having the same 
center line as the Impeller produces a parabola on the developed 
cylinder surface. Because of adherence to radial blade elements, 
the blade angle 8 increases from hub to shroud. 

The operating conditions for the Impeller were as follows: 


Equivalent impeller tip speed, Tj/ffd, ft/sec 1331 

Batio of specific heats, 7 ....... 1.4= 

Ererotatian at inlet ....■ 0 

Impeller equivalent weight flow, W A/fT/B, lb/sec 8.743 


Procedures for Numerical Solution 

The meridional plane of the impeller together with the estimated 
streamlines and resulting normals is shown in figure 1(a). Estimated 
streamlines are drawn in between the hub and the shroud; equal mass 
flow is assumed between the annular stream tubes. Lines normal to 
the streamlines are drawn in along the flow path. 

Radius of curvature . - The radius of curvature of the stream- 
line at any point in the meridional plane r c is obtained by 

r [l + (r- )!] 
r p 

»» j.n 

The slope of each streamline r* at each normal n is plotted 
against a xial distance in order to obtain the change in slope r” . 
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Numerical Integration * - The velocity q is found by numerically 
integrating equation (9) 


r 


a dn 


rn 


q. = e 


- 


uo 


-r 


a dn 


he 



( 9 ) 


The velocity is the estimated velocity at the huh* 

By plotting a against n (fig* 3(a)) 'where n is the distance 
along the normal from huh to shroud, the integral a dn can he 


evaluated hy the trapezoidal rule. A plot of h e 


-JC 


a dn 


against n can then he obtained (fig. 3(h)). and evaluated in a 
similar manner. As previously stated, the term du/dm in para- 
meter h is neglected in the first solution. The evaluations of 
these integrals allow the determination of the velocity q at 
any point n. 

The weight flow at any normal may he obtained from equation (10) 




2 * 


. yJO 


rfpgq cob P dn 


( 10 ) 


When a plot of rfpgq cos P against n 1 b used' (fig* 3(c)), the 


integral 


I 


rfpgq cos P dn can he evaluated hy the trapezoidal 


rule. If the weight flow does not equal the designated weight flow 
for the 'passage, a new value of the velocity q 0 is assumed and the 

solution is repeated. 


In order to satisfy continuity for each stream tube, new 
streamlines are drawn based on the velocities of the preceding 
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solution. Velocity q (fig* 4(a)) and density Pg (fig. 4(b)) 
are plotted against n. The radii at which the new streamlines 
cross the normals are obtained as follows: Erorn. 


and 


An 


r n+l ~ r n 
cos a 


/ 


Aw = 2itrfpgq. cos 6 An 


n+1 



Aw cos a 
Jtfpgq cos 0 


where cos a, f , p, q, and cos 0, are average values between 
streamlines based on estimated values of r n+1 . After the new stream- 
lines are determined, the normals are drawn. The graphical solution 
is then repeated for succeeding sets of new streamlines until con- 
tinuity is satisfied. 


O 


Accuracy 

The accuracy of a solution depends to a large extent upon the 
accuracy with which the streamline curvature Ib determined. Measure- 
ment of other geometrical quantities, of course, alters into the 
accuracy, but their measurement Is more precise. 

The sudden decrease in flow area (as much as 25 percent at the 
hub) from a point Just ahead of the blades at the inlet to a point 
Just inside the blades distorts the streamline picture. A tedious 
and detailed solution is necessary for a true picture of the dis- 
tortion in this region. Also, the assumption that the flow angle 
is the same as that of the blades at the inlet is questionable. 

A lengthy analysis at the inlet was therefore considered 
unjustifiable . 

The cross-sectional, area for the clearance is approximately 
4 percent near the inlet and affects the maximum flow. In the 
solution, the clearance area was treated as though the blades extended 
to the shroud because of continuity considerations. 
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Besults and Discussion 

O 


A solution of the flov conditions throughout the passage from 
hub to shroud Is shown in figure 5* 

Streamlines * - The flow streamlines in the meridional plane' 
are shown in figure 5(a). Although the solution has been made for 
equal quantities of flow between streamlines, little information 
can be derived from the streamline configuration, as shown, other 
than the direction of flow because the variation in distance 
between streamlines is due to a variation in blade thickness and 
radius as well as in velocity. Hear the outlet where blade thick- 
ness and radius have little effect, however, the streamline spacing 
indicates considerable variation in velocity across the passage. 

Telocity ratio . - Lines of constant velocity ratio Q relative 
to the impeller are shown on figure 5(b). Along the shroud (fig. 6), 
the velocity is near sonic at the inlet (Mach number, 0.973) with 
immediate deceleration just inside the impeller and continued 
deceleration to the outlet, a condition that is highly conducive to 
boundary-layer buildup and flow separation. A similar condition 
exists along the straight part of the hub (A to B, fig. 5(b)). 

Along the curved part of the hub (B to C), the velocity variation 
is small with a slight acceleration existing near the impeller outlet. 

The relative -velocity distribution aB obtained in the solution 
represents the mean velocity from blade to blade. According to the 
analysis of reference 1, the velocity along the suction surface is 
much larger than that along the pressure surface; the velocity varies 
Bomewhat uniformly across the passage. If the blade unloads at the 
outlet, as it must to satisfy the Kutta condition of tangency at the 
tip, the decelerations along .the suction surface of this blade are 
of necessity much greater than that shown far the average velocity 
(fig. 5(b)). Such adverse velocity gradients could easily cause 
considerable separation along the flow path. At the impeller outlet, 
these flow separations would result in poor flow conditions with 
consequent losses at the diffuser inlet. 

Pressure ratio . - lines of constant pressure ratio are shown 
on figure 5(c). The pressure lines approximately parallel the normal 
lines and indicate that static-pressure readings on the stationary 
shroud of this impeller give representative values for the static 
pressure across the passage from hub to shroud. In figure 7, 
experimental values of pressure ratio along the shroud are compared 
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with the theoretical values obtained in the numerical solution. 
Pressure ratios are shown for the theoretical weight flow and ‘brack- 
eting experimental weight flows. The experimental curve for the 
higher weight flow represents a point very near the maximum flow 
or flow cut-off line and shows a large pressure loss near the inlet. 
The experimental curve for the- lower weight flow represents an 
operating point near maximum efficiency and consequently approaches 
the theoretical conditions. No pressure drop occurs at the inlet, 
which is in agreement with the theoretical curve; however, the 
theoretical pressure increases more rapidly along the shroud because 
of inefficiency due to friction and flow separation, velocity 
Increase resulting from flow-area reduction due to boundary layer, 
and slightly higher theoretical tip speed. 

Maximum flow . - By making a series of estimates of velocity 
at the hub and obtaining the corresponding weight flows, the maximum 
theoretical weight flow across any normal can be determined. Focr 
the numerical example, the greatest restriction of flow occurred 
just Inside the inlet. The maximum flow over a range of tip speeds 
was determined by the consideration of this flow restriction. Inas- 
much as the radius of curvature was Infinite for the streamlines 
near the inlet in the complete numerical solution, it was taken as 
infinite in computing maximum weight flow at all speeds. A com- 
parison of the theoretical and experimental maximum flow over the 
speed range is shown In figure 8(a). The slightly higher theoretical 
values are in good agreement with the experimental values over a 
range of tip speeds up to 1400 feet per second (deviation of the 
order of 3 percent). At speeds above 1400 feet per second, the 
deviation increased to a maximum value of 7 percent at a tip speed 
of 1600 feet per second. 

The high relative velocities within the impeller, which caused 
choking, resulted from the large blade angle P near the Impeller 
Inlet. A large Inlet blade angle was required for shockless inlet 
conditions resulting from the consideration of the impeller mass 
flow and the rotational speed. - These high relative velocities, 
however, can be reduced by an appropriate selection of blade shape 
to give a more rapid decrease in £ near the impeller inlet. 

A comparison of P for the example impeller A and a second impel- 
ler B (reference 3) is shown in figure 8(b); P for each impeller 
is plotted against the distance ratio L along the shroud. Impel- 
ler B has the same hub -to - shroud profile and inlet blade angle p 
as 'that of impeller A; however, as indicated in figure 8(b), the rate 
of decrease in blade angle for impeller B is much larger near the 
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impeller inlet* The resulting higher flow capacity of impeller B 
(fig* 8(a)) may he attributed to the larger rate of decrease of 
blade angle and consequently lower relative velocities in the 
impeller passage near the inlet. 


STM4AEY OF RESULTS 

A method has been developed far a compressible-flow analysis 
between the hub and the shroud of mixed -flew impellers of arbitrary 
design* Axial symmetry was assumed, but the forces in the meridional 
plane that are derived from tangential pressure gradients were taken 
into account. Application of the method to a mixed -flow impeller 
for which experimental data are available yielded the following 
results : 

1* A large variation in velocity existed across the passage 
from hub to shroud. 

2. Adverse velocity gradients occurred along both hub and shroud 
surfaces. 

. 3. The theoretical pressure ratio increased more rapidly than 
the experimental along the shroud from inlet to outlet. Inefficiency 
of the actual compression process due to flow separation accounts 
for the greater part of this deviation. 

4. The maximum flow over a range of speeds as theoretically 
determined closely agreed with the experimental. 

5. For the example calculation, values of static pressure at 
any point along the shroud were representative of those across the 
flow path from hub to shroud. 


CONCLUSION 

Analysis of flow in the meridional plane of a mixed-flow impeller 
showed that the rotational farces, the blade curvature, and the hub- 
shroud profile can introduce adverse velocity gradients! along the hub 
and shroud surfaces. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, April 13, 1950. 
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APPENDIX A 
SYMBOLS 

The following symbols are used In the calculations and the 
figures : 

, cos*' 6 sin*' 8 dr sin*^ 8 . 

a parameter equal to c z sin a tan a 

r c r dn r 

b parameter equal to 2 oj sin 8 ^ + sin 8 tan a (zu) sin a + jj~ ^ 

r 

c speed of sound for stagnation conditions; ft/sec 
F force per unit mass. 

f factor accounting for area taken up by blade cross section 

g acceleration due to gravity, ft/sec 2 
H absolute total enthalpy, Btu/lb 

h static enthalpy, Btu/lb 

L ratio of distance m at any point to m at impeller exit 

m distance along streamline from inlet, ft 

n distance along normal to streamline in meridional plane, ft 
p static pressure, lb/sq ft. 

Q velocity ratio, q/cj. 

q velocity relative to impeller (fig. 1(b)), ft/sec 

r radius of particle from axis of impeller rotation, ft 

r radius of curvature of streamlines in meridional plane (positive 
c ■ in fig. 2(a)), ft 

r* slope of streamline in meridional plane .with respect to axis 
of rotation 
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r" change in slope r‘ with respect to axis of rotation 
t time, sec 

U actual impeller tip speed, ft/sec 

u. tangential velocity relative to impeller (positive in direction 
of rotation, fig. 2(c)), ft/sec 

v through-flow component of velocity (fig. 2(c)), ft/sec 

W total compressor flow rate, Xb/sec 

w flow rate, Ib/sec 

z distance in axial direction, ft 

■ 

a angle between tangent to streamline in meridional plane and 
axis of rotation, deg 

6 blade angle (fig. 1, negative for backward - swept blades), deg 

7 ratio of specific heats 

5 ratio of total pressure at inlet to standard sea-level pressure 

0 ratio of inlet stagnation temperature to standard sea -level 

temperature 

A. prerotation term equal to r^ua^ + u^) 

2 / 4 

P mass density, lb-sec /ft 

2 / 4 

p . total mass density, lb-sec /ft 

9 , blade angle in axial -tangential direction (negative for backward- 
swept blades, fig. 2(c)), deg 

w angular velocity of impeller, radians /sec 
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Subscripts : 

1 inlet conditions 

n indicates position along normal from hub to shroud or normal 
component 

r radial component 

z axial component 

v 

9 tangential component 

0,1 points at succeeding streamlines along normal 
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APPENDIX B 


DEVELOPMENT OF FORCE EQUATION IN HOB-TO -SHROUD PLANE FOR 

FLOW WHEN AXIAL SYMMETRY IS ASSUMED 

Lorenz's equations of motion relative to the rotating impeller, 
in which axial symmetry is assumed (reference 3), are 


F. 


1 dp . d^ (wr + u)‘ 

"* T _ ** 11 


r P ^ dt 2 


F _ i & + 

2 p at 2 


F du u dr _ dr 

3 _ + _ _ + 2 ( 1 ) — 

9 dt r dt dt 

where the "blade force F is the force per unit mass. 

From figure 2(a), 

F = F sin a - F cos a 
■ n z r 

where F may he evaluated from F_ (appendix C) 
n , o 


F. 


1 dp . d 2 z . 

= — ~ sin a + sin 


Q P dz 


1 dp d^r 

a tt 5 - cos a 


dt 4 



dr 

dt 


dz 


dt 

d^r 

dt 2 

= V cos 


p dr 


= v sin a 


= v cos a 


da 

dt 


dt £ 


(car+u) 2 

cos a + — cos a 


(HI) 


dv 

dt 


(Bla) 
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d^z da dv 

=-z = - v sin otTr + cos a rr 

at 2 4t dt 

& = &«£ + &£ = .& sin o. + |E cos a 
dn dz dn or dn oz or 

Substituting equations (Bla) to (Blc) in equation (Bl) gives 


F « . iS£ - y + fegl 

*n p dn dt r 


cos a 


Item figure 2(a), 


da = 


v dt 


Equation (B2) therefore reduces to 


if . (tur+tt) 2 dr g 

r r dn ” J n 

c 


(Bib) 

(Blc) 

(B2) 


1 dn 
p dn 


+ 


(B3) 
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APPENDIX C 

PROCEDURE FOR EVALUATING BLADE FORCE F FOR 

n 

IMPELLER OF ARBITRARY BIADE SHAPE HAYING 
RADIAL BIADE ELEMENTS 

For an impeller with, radial "blade elements, a true picture of 
the "blade curvature lies on a circular cylinder having the same 
center line as the impeller. The "blade force F normal to the "blade 
(fig. 2(h)) has the component F 0 , 'which may he determined from time 

rate of change of angular momentum in the direction of rotation. The 
component F 0 as given in appendix B mates the angle q> with the 
force F. 

From figure 2(c), 


and 


F = -F- tan cp sin a 
n , y 


tan 


<P = 


u 

v cos a 


tan p 
cob a 


o 


giving 


From appendix B, 


F = -F Q tan 6 tan a 
n (7 


du u dr „ dr 

F_ = Tr + -vr + 2U)=~ 
0 dt r dt dt 


Substituting dm = v dt and dr/dm = sin a gives 
Fq «= v ~ sin a + 2 U) sin 

F = -v tan B tan a (~ + - sin a + 2co sin a) 
n \am r / 
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3P 

n 


a ^2 


-v tan 3 tan a ( 20) sin a + 


du\ _ uv 
dm j ~ r 


sin a tan a tan 3 


From figure 1, 


v =■ q cos 3 


u = q sin 3 


F„ = -q sin 3 tan a 
n 


2u) Bin a. + 


du\ 

dm/ 



P 


sin a tan a 
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Meridional 

streamline 




(a) Meridional plane components. 






(!) Components of force normal 
to Hade. 


(o) Orientation of tp with 

respect to p. (Bashed 
lines lie in meridional 
plane) 


Figure 2. - Quantities used in development of farce equations. 



Distance along normal to streamline, n, ft 


(o) Weight -flour factor. 


(a) Parameter a. ("b) Exponential h e (o) Weight- 

Plgure 3. - Plots use! for numerical integration of equations (10 ) and (11). 


aui 






Distance along normal to streamline, 
n, ft 

(t) Density. 

Figure 4. - Plots used for computing 
new streamlines. 







1336 




HACA nr 2165 


Hatlo of static pressure 
to Inlet total pressure 
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Figure 6. - Comparison of relative velocity and "blade 
angle along shroud.. 




Ratio of static pressure t6 inlet total pressure 
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Maximum equivalent weight flow. 

Blade angle, P, deg lb/sec 
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600 800 1000 1200 1400 1600 

Equivalent tip speed, Tj//<f, ft/Bec 


(a) Comparison of maximum flow capacity . 



0 .2 .4 .6 .8 1.0 

Distance ratio along shroud, L 


(b) Comparison of blade-angle decrease from inlet 
to outlet at shroud. 

Figure 8. - Comparison of change in blade angle ma-Ti-m-nm 
flow capacity of two impellers having same hub -shroud 
profiles and blade inlet angles. 
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